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Abstract 

We introduce a class of Jacobi operators with discrete spectra which is character- 
ized by a simple convergence condition. With any operator J from this class we 
associate a characteristic function as an analytic function on a suitable domain, 
and show that its zero set actually coincides with the set of eigenvalues of J in 
that domain. Further we derive sufficient conditions under which the spectrum 
of J is approximated by spectra of truncated finite-dimensional Jacobi matrices. 
As an application we construct several examples of Jacobi matrices for which the 
characteristic function can be expressed in terms of special functions. In more 
detail we study the example where the diagonal sequence of J is linear while the 
neighboring parallels to the diagonal are constant. 
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1 Introduction 

In this paper we introduce and study a class of infinite symmetric but in general 
complex Jacobi matrices J characterized by a simple convergence condition. This 
class is also distinguished by the discrete character of spectra of the corresponding 
Jacobi operators. Doing so we extend and generalize an approach to Jacobi matrices 
which was originally initiated, under much more restricted circumstances, in [20J. We 
refer to [TH] for a rather general analysis of how the character of spectrum of a Jacobi 
operator may depend on the asymptotic behavior of weights. 

For a given Jacobi matrix J7", one constructs a characteristic function Fj(z) as an 
analytic function on the domain Cq obtained by excluding from the complex plane 
the closure of the range of the diagonal sequence A of J . Under some comparatively 
simple additional assumptions, like requiring the real part of A to be semibounded or 
J to be real, one can show that J determines a unique closed Jacobi operator J on 
£ 2 (N). Moreover, the spectrum of J in the domain Cq is discrete and coincides with 
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the zero set of Fj(z). When establishing this relationship one may also treat the poles 
of Fj(z) which occur at the points from the range of the sequence A not belonging 
to the set of accumulation points, however. In addition, as an important step of the 
proof, one makes use of an explicit formula for the Green function associated with J. 

Apart of the localization of the spectrum we address too the question of approx- 
imation of the spectrum by spectra of truncated finite-dimensional Jacobi matrices. 
For bounded Hermitian Jacobi operators the problem has been studied, for example 
in [TTj [3J [13]. We are aware of just a few papers, however, bringing some results in this 
respect also about unbounded Jacobi operators [T^ [T2~| ITS] . Our approach based on 
employing the characteristic function makes it possible to derive sufficient conditions 
under which such an approximation can be verified. This result partially reproduces 
and overlaps with some theorems from |12j . 

The characteristic function as well as numerous formulas throughout the paper are 
expressed in terms of a function, called 5, defined on a subset of the space of complex 
sequences. In the introductory part we recall from [20] the definition of J and its basic 
properties which are then completed by various additional facts. On the other hand, we 
conclude the paper with some applications of the derived results. We present several 
examples of Jacobi matrices for which the characteristic function can be expressed in 
terms of special functions (the Bessel functions or the basic hyper geometric series). A 
particular attention is paid to the example where the diagonal sequence A is linear while 
the neighboring parallels to the diagonal are constant. In this case the characteristic 
equation in the variable z reads J„ 2 (2u>) = 0, with w being a parameter, and our main 
concern is how the spectrum of the Jacobi operator depends on w. 



2 The function $ 

2.1 Definition and basic properties 

Let us recall from [2D] some basic definitions and properties concerning a function $ 
defined on a subset of the linear space formed by all complex sequences x = {xk}kLi- 
Moreover, we complete this brief overview by a few additional facts. 

Definition 1 . Define $ : D — > C, 

oo oo oo oo 

#0*0 = 1 + 5^(-l) m Yl Yl Y X kl X kl+1 X k2 X k2+1 . . -X km X km+1 , (1) 

m=l fci = l A)2=fci+2 fc m =fc m _i+2 

where 

{oo 
{x k }lZi CC;Y \ x kXk+i\ < oo 
k=l 

For a finite number of complex variables we identify ^{x\, X2, • • • , x n ) with $(x) where 
x = (xi, X2, . . • , x n , 0, 0, 0, ... ). By convention, we also put 5(0) = 1 where is the 
empty sequence. 
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Let us remark that the value of $ on a finite complex sequence can be expressed 
as the determinant of a finite Jacobi matrix. Using some basic linear algebra it is easy 
to show that, for n G N and {xj}™ =1 C C, one has 



where 



$(xi,x 2 , ...,x n ) = detX n 

/l 21 

x 2 1 x 2 



(2) 



X„ 



\ 



X n —\ 1 X n —\ 
\ x n 1 J 

Note that the domain D is not a linear space. One has, however, £ 2 (N) C D. In 
fact, the absolute value of the mth summand on the RHS of ([I]) is majorized by the 
expression 



^ , i 1 (sr 

;eN m ' \j=l 



XjXj + \ I 



keN 
k\<k2<---<k 



Hence for x G D one has the estimate 

|5(x)|<exp(]T 



X k X k+ \\ 



Furthermore, $ satisfies the relation 
$(x) = $(xx, ...,x k ) $(T k x) - !$(xi, x k -i)x k x k+1 $(T k+1 x), k = 1, 2, 



(3) 



(4) 



where x G D and T denotes the truncation operator from the left defined on the space 
of all sequences, T({x n }^ =1 ) = {z n+ i}™ =1 . In particular, for k — 1 one gets the rule 

$(x) = 3{Tx) — x 1 x 2 $(T 2 x). (5) 

In addition, one has the symmetry property 

X 2 , ■ ■ ■ , Xk-l, X k ) = $(Xk, Xk-l, ■ ■ ■ , X 2 , Xi). 

If combined with (J3J), one gets 

$(xi,x 2 , x k+ i) = 3(xi, x 2 , . . . , x k ) - x k x k+1 $(xx,x 2 , x k -i). (6) 
Lemma 2. For x G D one has 

lim $(T n x) = 1 

n— >oo 

and 



lim x 2 , ...,x n ) = $(x). 



(7) 
(8) 
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Proof. First, similarly as in ([3]), one gets the estimate 

/ oo 

\${T n x) - 1| < exp J2 \ x kX k+1 \ ) - 1. 



v fc=n+l 



This shows 0. 

Second, in view of 01]), the difference \$(x) — $(xi,X2, ■ ■ ■ ,x n )\ can be majorized 
by the expression 



oo 

x k x k+l | 



|1 -£(T n x)|exp| ^\x k x k+l \ J + |x n x n+ i|exp| 2^ 

\k=l J \ fc=l 

From here one derives the (rather rough) estimate 

(oo \ oo 

x k x k+l 
| ^ \ X k%k+l\- (9) 
k=l J k=n 

This shows ©. □ 

Proposition 3. The function $ is continuous on £ 2 (N). 

Proof. If x G £ 2 (N) C -D then from ([5]) one derives that, for any nGN, 



\$(x)-!$(x 1 ,x 27 ...,x n )\ < 2 exp(2||x|| 2 ) ||(/-P n -i; 



x\\ 2 



where P m stands for the orthogonal projection on £ 2 (N) onto the subspace spanned by 
the first m vectors of the standard basis. From this estimate and from the fact that 
$(xi, X2, ■ ■ ■ , x n ) is a polynomial function the proposition readily follows. □ 



2.2 Jacobi matrices 



Let us denote by J an infinite Jacobi matrix of the form 

[XiWi \ 

J = 



Wi 

V-2 



W 2 
A3 



w 3 



V 



•7 



where w = {w n } 



n=li 



{v n }™ = i C C \ {0} and A = {A„}^° =1 C C. Provided any of 
the sequences is unbounded it is reasonable to distinguish in the notation between J 
and an operator represented by this matrix. Such an operator J need not be unique, 
as discussed in Subsection 13.21 Further, by J n we denote the nth truncation of J~, i.e. 



/ Ai Wi 

vi A 2 w 2 



J,, 



\ 



V n - 2 A n _i W n -i 

v n -i A n 



(10) 
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As is well known and in fact quite obvious, any solution {x k } of the formal eigen- 
value equation 

Ai^a + Wtx 2 = zx u v^x^i + X k x k + w k x k+ i = zx k for k > 2, (11) 

with z G C, is unambiguously determined by its first component X\. Consequently, 
any operator J whose matrix equals J may have only simple eigenvalues. 

We wish to show that the characteristic function of a finite Jacobi matrix J n can be 
expressed in terms of J- To this end, let us introduce the sequences {7^}/! =1 defined 
recursively by 

it = 1, 7^+i = w k /^ and 7^ +1 = v k /-y£, k > 1. (12) 
More explicitly, the sequence {7/7}fc =1 can be expressed as 

As for the sequence {7^)^=1, the corresponding expressions are of the same form 
but with w being replaced by v and vice versa. Note that if v k = w k for all k = 
1, 2, . . . , n — 1, then 7^ = j k for all = 1, 2, . . . , n. 

Proposition 4. Let {7^}^=! be the sequences defined in HM) . Then the equality 

d*w. - = (n<A, (13) 

/ioWs /or all z & C (after obvious cancellations, the RHS is well defined even for 
z = X k ). 

Proof. Put Afc = A / t/7 A T7 A l". As remarked in [2"U| Remark 24], the Jacobi matrix J n 
can be decomposed into the product J n = G~J n G+ where = diag(7^, 7^, . . . , 7^) 
are diagonal matrices, and J n is a Jacobi matrix whose diagonal equals the sequence 
(Ai, A2, • • • , A n ) and which has all units on the neighboring parallels to the diagonal. 
The proposition now readily follows from this decomposition combined with ([2]). □ 

Moreover, with the aid of (lT5j) and using some basic calculus from linear algebra 
one can derive the following formula for the resolvent. 

Proposition 5. The matrix entries of the resolvent R n (z) = (J n — zJ n ) -1 , with z G 
C \ spec(J n ) ; may be expressed as (1 <i,j <n) 



(max(i.j') \ 
n (* - a <) ( m ) 
l=mm(i,j) J 



A/ — Z I j -t j \|A/ — z \ , rn r • - \ 1 -1 / VIA; — Z 1 ; 1 

1 s 1=1 I \ v. 1 J i=max(tj)+l / v v. t j 1=1 



where 



lU^w h ifi<j, 
= { 1, ifi = j, 

Y[]l)vu %fi>j. 



In the remainder of the paper we concentrate, however, on symmetric Jacobi ma- 
trices with v = w, i.e. we put 



J 



( Ai wi \ 

w 1 \ 2 w 2 

w 2 A 3 w 3 

V '■) 



where A = {X n }^i C C and w = {w n }^ =1 C C \ {0}. In that case some definitions 
introduced above simplify. First of all, one has 7^ = 7^ = 7& where 



k-l 



72fc-l 



n 



W 2 j 



1 W 2j-l 



72/b = wi 



w 2j+1 



k-l 



k = 1,2,3, 



Then 7 fc 7fc+i = w k . 



2.3 More on the function $ 

In [20] one can find two examples of special functions expressed in terms of The 
first example is concerned with the Bessel functions of the first kind. In more detail, 
for w, v G C, v ^ — N, one has 



w 



r(i/ + l) 

Notice that jointly with ([3]) this implies 



w 



u + k 



\JJ2w)\ < 



w 



vfc = l 



k=l 



W 



(15) 



r(i/ + i) 

In the second example one shows that the formula 

^m(2m-l) r np' m 



{u + k){u + k + 1) 



(16) 



3 r «**-ML) = i+E(- 1 ) m 7i 



rn=l 



(l-t 2 )(l-t 4 )...(l-t 2m ) 



o&GO;* 2 ,-^ 2 ) (17) 



holds for i, w G C, |t| < 1. Here o0i is the basic hypergeometric series (also called 
q-hypergeometric series) being defined by 

o0i(; 6; g, z) = ^ " — — — ~ zfc 



fc=0 
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{q;q)k{b; q) k 



and 

fe-i 

(a;q) k = Y[(l-aq j ), k = 0,1,2, ... , 

3=0 

is the g-Pochhammer symbol; see jH] for more details. 

In this connection let us recall one more identity proved in fZU\ Lemma 9] , namely 

Ui $ (U2, U 3 ,..., U n ) # (V l, V 2 ,..., V n ) -V!$ (Ul, U 2 ,..., U n ) # («2, ^3, • • • , V n ) 

= ^2 [ II UkVk ) ( u i ~ v ^ $ ( u j+ii M i+ 2 ' • • • ' Un ) $ ( v i+i' v i+^ ■■■i v n) ■ 
j=i \k=i J 

For the particular choice 

w w 
u k = — —r , V k = — — r , 1 < k < n, 
fj, + k v + k 

one can consider the limit n — > oo. Using (Il5|) and ( ITS]) one arrives at the equation 

oo 

J M (2w) J, + i(2w) - J M+ i(2«;) J„(2tw) = ^ J M+i (2«;) J, +J (2w). (18) 

i=i 

Definition (JT]) can naturally be extended to more general ranges of indices. For 
any sequence {x n }%L Nl , N 1} iV 2 G Z U {-oo, +00}, N x < N 2 + 1, (if Nx = N 2 + 1 G Z 
then the sequence is considered as empty) such that 

N2-1 

\xkXk+i\ < 00 

k=N 1 

one can define 

00 

d{[xk}kl Nl ) = 1 + ^h 1 )™ x kl Xk 1+ ix k2 Xk 2+ i . . . x km x km+ i 

m=l kel(N 1 ,N 2 ,m) 

where 

X(Nx, N 2 , m) = {ke Z m ; kj + 2 < k j+1 for 1 < j < m - 1, Ni < h, k m < N 2 } . 
With this definition one can generalize the rule (j2J). Now one has 

^({a*}^) = S( WLjvJ s({^}f= n +i) - x„x n+ i5({x4;4) s({^}f=n+ 2 ) 

(19) 

provided n G Z satisfies iVi < n < N 2 . 

This extension also opens the way for applications of the function J to bilateral 
difference equations. Suppose that sequences {w n }'^L_ oo and {Cn}^L_oo are such that 
w n 7^ 0, C, n 7^ for all n and 
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E 



k=— oo 



< OO. 



CfcCfc+1 

Consider the difference equation 

W n U n+ i - CnU n + %-l«n-l = 0, 71 E Z. (20) 

Define the sequence {V n } neZ by V = 1 and P n+ i = (w n /( n+ i)V n for all n. Hence 

o 



for n > 0, P = 1, P n 



n 



c* 



Wk-1 



for n < 0. 



fc=l " fc=n+l 

The sequence {7 n } ne z is again defined so that 71 = 1 and 7„7 n +i = w n for all iiGZ, 
Hence 



fe-i 



72fc-l 



n 



w 2i 



1 W 2j-1 



, 72fc = Wi 



fc-1 



, for = 1, 2, 3, ... , 



and 





T2k-i = n 

i=fc 



W 2j-1 

w 2j 



72fc = 



n 



w 2i 



Then the sequences {f n } n& and {fi'„} neZ , 



for fc = 0,-1,-2, 



fn = V n d 



2 \ 00 



fc=n+l 



1 



9r, 



W n -{P, 



n-l 



1 ~\ n—1 
ill 
Ck J k= _, 



(21) 



represent two solutions of the bilateral difference equation (120]) . 

For two solutions u = {u n } nl - z and v = {v n } ne z °f (EHJ) the Wronskian is introduced 

as 

W{u, v) = w n {u n v n+ i - U n+1 V n ) . 

As is well known, this is a constant independent of the index n. Moreover, two 
solutions are linearly dependent iff their Wronskian vanishes. For the solutions / and 
g given in (|2ip one can use f lT9|) to evaluate their Wronskian getting 



W,</) = 3 



9 -\ 00 



One may also consider an application of a discrete analogue of Green's formula to 
the solutions (T2"Tj) [2J. In general, suppose that sequences {w n }^lo and {f n }^L solve 
respectively the difference equations 

W n U n+ i - Ci 1} «n + W„-lW„-l = 0, W n f n+ 1 - Cn ] V n + W n -iV n -i = 0, U G N. (22) 

In that case it is well known and easy to check that 

n 

^2 (ti^ ~ U 3 V 3 = W ° ( U ° Vl ~ UlV °^ ~ Wn ( U ™ V n+l ~ Un+lVn) ■ (23) 

i=i 
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Proposition 6. Suppose that the convergence condition 



fc=i 



CfcCfcH 



< OO 



is satisfied for the both difference equations in ( d|| ). Moreover, assume that 



sup 

n>l 



w; 



a(1)a(2) 
S>« Sn+l 



< oo and sup 

n>l 



a(2)a(1) 



+ 1 



< OO. 



Then the corresponding solutions f^\ f^ from ![21\) fulfill 

oo 

£ (Cj 1} - Cf ) if if = -o (/ (1) A (2) - /fif 



Proof. In view of (123 j) it suffices to show that 

hm «; B /W/Si = lim ^tfi/f = 0. 

ra— >oo n— >oo 

By the convergence assumption, for all n > one has 



ki< Wic^iicSii, w<^vicrnc 



(2) 1 1/-( 2 ) I 
n+ll 



Using ([2]), after some straightforward manipulations one gets the estimate 



w f (1) f (2) 

w nJ n Jn+l 



exp 



/ oo 

\fc=l 


w k 


+ 




\ "0 

/ fc=l 


">k-l 






a(2)a(2) 


a(1)a(2) 



i r (l) r (2)|l/2 



aWa(2) 
Sn+1 



1/2 



This implies 



(24) 



(25) 



□ 



In the literature on Jacobi matrices one encounters a construction of an infinite 
matrix associated with the bilateral difference equation (!2Up |21[ § 1.1], |1U[ Theo- 
rem 1.2]. Let us define the matrix 3 with entries 3(m,n), m,n £ Z, so that for every 
fixed m, the sequence u n = 3(m, n), n £ Z, solves ( 120]) with the initial conditions 
3(m, m) = 0, 3(m, m + 1) = l/iu m . 

Using ([6]) one verifies that, for m < n, 



J(m,n) 



; ( H 



j=m+l J , 



7m+l 7m+2 7n-l 
Cm+1 Cm+2 Cn-1 



Moreover, it is quite obvious that, for all m,n £ Z, 

1 



3(m, n) 



W(«,v) 



[u m v n v m u n ) , 
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where {u n }, {v n } is any couple of independent solutions of ([20]) . Hence the matrix Z 
is antisymmetric. It also follows that, Vm,n,k,£ G Z, 

3(m, k)3(n, £) - 3(m, £)3(n, k) = 3(m, n)Z(k, £). 

Example 7. As an example let us again have a look at the particular case where 
w n = w, ( n = v + n for all n G Z and some w^eQw/O^^Z. One finds, with 
the aid of f lT5|) . that the solutions f l2T]) now read 

(-l) n 7T 

/ n = T(z/ + 1) w v J u+n (2w), g n = — — — — — w v J- v - n (2w). 

sm(7ri/)l (z/ + 1) 

Hence the Wronskian equals 

7TW ( { W 1 °° \ 

>V /,<? = -7-^ (—J u (2w)J- u -i(2w) - J v+1 (2w)J„ v (2w)) = d( — — 
Recalling once more ( 115]) we note that the RHS equals 



lim $1 < } = lim > i — / — i >—w 2n = 1. 

n^oo \[p - N + nj n= J n-^oo n! r(z/-iV + n+l) 

Thus one gets the well known relation [I] Eq. 9.1.15] 

J„ +l (2w) J_ u (2w) + J v {2w)J^{2w) = - S -^™1 . (26) 

Concerning the matrix 3, this particular choice brings us to the case discussed in 
[20"] Proposition 22]. Then the Bessel functions Y n+U (2w) and J n+U (2w), depending 
on the index nGZ, represent other two linearly independent solutions of (1201) . Since 
dJ Eq. 9.1.16] 

J v+1 {z)Y v {z) - J u (z)Y u+1 (z) = — 

TTZ 

one finds that 

3(m,n) = tc (Y m+u (2w)J n+v (2w) - J m+u (2w)Y n+v (2w)) . 
Moreover, for a = m + \i and k = n — m > one has 



r(er + k) I f w 



J a+k (2w)Y a (2w) - J a (2w)Y a+k (2w) = ^— 5 



fc-l x 



7TU7 fc r((J + 1) V V CT + J J J = 1 

- and w n = w, Wn G N, in equs 
one verifies that (|24|) holds true and reveals this way once more the identity (1 18)) 



Finally, putting Cn = /U + n, ^ 2 ' ) = i/ + n and u>„ = w, Vn G N, in equation 
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3 A class of Jacobi operators with point spectra 
3.1 The characteristic function 

Being inspired by Proposition H] and notably by equation ()13j) . we introduce the (renor- 
malized) characteristic function associated with a Jacobi matrix J as 

™==<{AD (27) 

It is treated as a complex function of a complex variable z and is well defined provided 
the sequence in the argument of $ belongs to the domain D. Let us show that this is 
guaranteed under the assumption that there exists zq G C such that 



n=l 



(A„ — 2 )(A n+ i — Zq) 



< oo. (28) 



For A = {Anj^^ let us denote 

C£ := C \ {A„; n G N}. 

Clearly, 

{A n ; jigN} = {A n ; n G N} U der(A) 

where der(A) stands for the set of all finite accumulation points of the sequence A (i.e., 
der(A) is equal to the set of limit points of all possible convergent subsequences of A). 

Lemma 8. Let condition [2$) be fulfilled for at least one zq G Cq. Then the series 

oo 2 

V ^ (29) 

l^{X n -z){\ n+x -z) 1 > 

converges absolutely and locally uniformly in z on Cq. Moreover, 

Mz G C* lim \ \ = Fj ( z ) : (30) 

and the convergence is locally uniform on Cq. Consequently, Fj(z) is a well defined 
analytic function on Cq. 

Proof. Let K C Cq be a compact subset. Then the ratio 

| A n — Zq\ ^ ^ \z — Zq\ 



| An Z | An ^| 

admits an upper bound, uniform in z G K and nGN. The uniform convergence on 
K of the series fl29l) thus becomes obvious. 
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The limit fl30l) follows from Lemma [2j Moreover, using f )30|) and also 



one 



has 



Ik 



A* 



oo 

^E 



n 

fc=l. 
2 



(A; - z)(\i+i - z) 



oo 

^E 

Z=n 
oo 



7fc 



exp 



E 

,fc=i 



A fc - z 



■w 



7k 



k=l, 



Afc 



fe=i. 



(A fc - z)(A fc +i - z) 



From this estimate and the locally uniform convergence of the series ( 129]) one deduces 
the locally uniform convergence of the sequence of functions (150]) . □ 

By a closer inspection one finds that, under the assumptions of Lemma [HI the 
function Fj(z) is meromorphic on C\der(A) with poles at the points z = \ n for some 
nGN (not belonging to der(A), however). For any such z, the order of the pole is less 
than or equal to r(z) where 

oo 
k=l 

is the number of members of the sequence A coinciding with z (hence r(z) = for 
z G Cq). To see this, suppose that r(z) > 1 and let M be the maximal index such 
that \m = z. Using (jlj) one derives that, for u G Cq, 



Fj(u) 




+ d 



. M \ 




I 


















"J 


r 

n= 





7; 



Ar. 



n=M+l 



7m7m+i 



z)(A 



A/+1 



5 



7« 



A, 



n=M+2 



The RHS clearly has a pole at the point u = z of order at most r(z). 



3.2 The Jacobi operator J 

Our goal is to investigate spectral properties of a closed operator J on £ 2 (N) whose 
matrix in the standard basis coincides with J . Provided the Jacobi matrix does 
not determine a bounded operator, however, there need not be a unique way how to 
introduce J. But among all admissible operators one may distinguish two particular 
cases which may respectively be regarded, in a natural way, as the minimal and the 
maximal operator with the required properties; see, for instance, [3]. 

Definition 9. The operator J max is defined so that 

Dom(J max ) = {y G £ 2 (N); Jy G £ 2 (N)}, 

and one sets J ma , x y = Jy, G Dom J max . Here and in what follows Jy is understood 
as the formal matrix product while treating y as a column vector. To define the 
operator J min one first introduces the operator J so that Dom(J) is the linear hull of 
the standard basis, and again Jy = Jy for all y G Dom(J). J is known to be closable 
[1], and J m in is defined as the closure of J. 
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One has the following relations between the operators J m i n , ./max and their adjoint 
operators [4, Lemma 2.1]. Let J H designates the Jacobi matrix obtained from J 
by taking the complex conjugate of each entry. Then J t * in = J 1 ^, J^ax = Jmm- m 
particular, the maximal operator J max is a closed extension of J mm . It is even true that 
any closed operator J whose domain contains the standard basis and whose matrix 
in this basis equals J fulfills J m j n C J C J max - Moreover, if J is Hermitian, i.e. 
J = J H (which means nothing but J is real), then J^ in = J max D J m i n . Hence J m i n 
is symmetric with the deficiency indices either (0, 0) or (1,1). 

We are primarily interested in the situation where J min = J max since then there 
exists a unique closed operator J defined by the Jacobi matrix J~, and it turns out that 
the spectrum of J is determined in a well defined sense by the characteristic function 
Fj(z). If this happens J is sometimes called proper [4]. 

Let us recall more details on this property. We remind the reader that the orthog- 
onal polynomials of the first kind, p n (z), are defined by the recurrence 

W„_lP n _l(z) + \ n Pn{z) + W n p n+1 (z) = Zp n (z), 71 = 1, 2, 3, ... , 

with the initial conditions p (z) = 1, pi(z) = {z—\i)/w\. The orthogonal polynomials 
of the second kind, q n (z), obey the same recurrence but the initial conditions are 
qo(z) = 0, qi(z) = l/u>i; see [21 E]. It is not difficult to verify that these polynomials 
are expressible in terms of the function $ as follows: 

and 

Wi \i=2 Wk J V U* - Z J 1=2/ 

The complex Jacobi matrix J is called determinate if at least one of the sequences 
p(0) = {p n (0)}^ = Q or q(0) = {g„(0)}^L is not an element of £ 2 (Z + ). For real Jacobi 
matrices there exits a parallel terminology. Instead of determinate one calls J limit 
point at +oo, and instead of indeterminate one calls J limit circle at +oo, see |2T| 
p. 48]. According to |2"2"| Theorem 22.1], J is indeterminate if both p(z) and q(z) are 
elements of £ 2 for at least one zgC, and in this case they are elements of £ 2 for all 
z G C. For a real Jacobi matrix J one can prove that it is proper if and only if it 
is determinate (or, in another terminology, limit point), see [21 pp. 138-141] or [2T| 
Lemma 2.16]. 

For complex Jacobi matrices one can also specify assumptions under which J min = 
J max . In what follows, p(A) designates the resolvent set of a closed operator A. Con- 
cerning the essential spectrum, one observes that spec ess ( J m in) = spec ess ( J max ) [U 
Eq. 2.10]. Hence if p(J ma x) ^ then spec ess ( J m i n ) ^ C. Moreover, in that case J 
is determinate [H Theorem 2.11 (a)] and proper [H Theorem 2.6 (a)]. This way one 
extracts from [4] the following result. 

Theorem 10. If p{J max ) ^ then Jmin = Jmax- 
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3.3 The spectrum and the zero set of the characteristic func- 
tion 

Let us define 

VJ) ■= [ z e £ \ der(A); lim (u - z) r{z) Fj{u) = o) . (31) 

Of course, 3(c7) H Cq is nothing but the set of zeros of Fj(z). Further, for k G Z + 
and z e C \ der(A) we put 



where one sets u>o = 1. One observes that for k > M, where M = M z is either the 
maximal index, if any, such that z = A a/, or M = otherwise, 

6W= n«*-(n(-*.)) HfeLJ- (33) 

Proposition 11. Let condition ( QSj) 6e fulfilled for at least one z E Cg. // 

fo(z) = lim (u - z) r <*>.F>(u) = 
/or some z G C \ der(A), £/ien z is an eigenvalue of J max and 

:=(eiW,6W,6W,...) 

is the corresponding eigenvector. 

Proof. Using (jSJ) one verifies that if £,o(z) = then the column vector £(z) solves the 
matrix equation J~£(z) = z£(z). To complete the proof one has to show that does 
not vanish and belongs to £ 2 (N). 

First, we claim that 7^ 0. Suppose, on the contrary, that £1(2:) = 0. Then 

the formal eigenvalue equation (which is a second order recurrence) implies £(z) = 0. 
From ([23]) it follows that 



for all k > M = M z . This equality is in contradiction with (171) , however. 

Second, suppose z ^ der(A) is fixed. By Lemma [SJ there exists N E N, N > M, 
such that 

K| < |A n - z||A n+ i - z\/2, Vn > N. 

Let us denote 

N N 



C = l[\w l _ 1 \ 2 } [ \z-X t \ 



1=1 1=1 
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Using also ((3]) one can estimate 



oo oo k k / f 2 } 00 \ 

Ei&wi^Eni^i'ni-^Kfc}, J 

k=N k=N 1=1 1=1 \K i J l=k+l/ 



1=1 



< Cexp|2 



k=N+l 



(Afe - z)(X k+1 - z) 



oo k 1 1 

n 

k=N l=N+l x 



z - X1-1 



z — Xi 



Since lAt- — z\ > r for all k > M and some r > 0, the RHS is finite. 



□ 



Further we wish to prove a statement converse to Proposition [TTJ Our approach is 
based on a formula for the Green function generalizing a similar result known for the 
finite-dimensional case; see (I14j) . 



Proposition 12. Let condition < fl?g)) 6e fulfilled for at least one z G Cq. If z E 

C\der(A) does noi belong to the zero set^(J') then z G p(J ma x) an d the Green function 
for the spectral parameter z, 

G(z; := (e h (J max - z^ej), i,j G N, 

(a matrix in the standard basis) is given by the formula 



G(z;i,j) 



1 



w 



max(ij') 

n 



max(i,i') \ , 

min(i 1 



z-X 



li 



X, 



1=1 




(34) 



i=max(i ,j)+l. 



If 



In particular, for the Weyl m-function one has 



m(z) := G(z;l,l) 



Ai — z 



li 



A; — Z 



1=2 



li 



Xi - z 



Xi-z 



oo \ — 1 



1=1 



1=1 



(35) 



If, in addition, \X n \ — > oo as n — > oo i/ien der(A) 
resolvent (J max — z)^ 1 is compact. 



and for every z G C\3(J7"), the 



Proof. Denote by R(z)ij the RHS of (134"]) . Thus R(z) is an infinite matrix provided its 
entries R(z) it j, i,j G N, make good sense. Suppose that a complex number z does not 
belong to 3{J) Uder(A). By LemmaEl in that case the RHS of ( 13~4"|) is well defined. By 
inspection of the expression one finds that this is so even if z happens to coincide with 
a member A^ of the sequence A not belonging to der(A), i.e. the seeming singularity 
at z = Afc is removable. For the sake of simplicity we assume in the remainder of 
the proof, however, that z does not belong to the range of the sequence A. The only 
purpose of this assumption is just to simplify the discussion and to avoid more complex 
expressions but otherwise it is not essential for the result. 
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First let us show that there exists a constant C, possibly depending on z but 
independent of the indices i, j, such that 



\R(z)i,j\ < C2- |i ~ J 'l, Vz,j G N. 



(36) 



To this end, denote 



t = inf{|z - A n |; n E N} > 0. 
Assuming (128]) . one can choose no G N so that, for all n > no, 

\w n \ 2 < |A„ - z\ |A n+ i - z|/4. 

Let us assume, for the sake of definiteness, that i < j. Again by 



(37) 



Ti 



A; - 2 



Z=l 



51 



If 



A/ - 2 J , 



=.7+1, 



7, 



and 



<Ci, 



for all i, j. It remains to estimate the expression 



| A, - z\ 



n 

l=i 



Wi 



Xi-z 



(38) 



We distinguish three cases. For the finite set of couples i, j, i < j < no, (13"8"]) is 
bounded from above by a constant C 2 . Using f l3T|) . if i < n < j then (I3"8]) is majorized 

by 



Co 



Ano % 



n 



A; - Z 



< Csr-^IA^-zl 1 /^-^" . 



Similarly, if n < % < j then (I3"5]) is majorized by r _1 2 _:;+i . From these partial upper 
bounds the estimate (136]) readily follows. 

From (136]) one deduces that the matrix R(z) represents a bounded operator on 
£ 2 (N). In fact, one can write R(z) as a countable sum, 



i2(z)=£i2(z; S ), 



(39) 



where the matrix elements of the summands are R(z; s)ij = R(z)ij if i — j = s and 
i?(z; s)jj = otherwise. Thus R(z; s) has nonvanishing elements on only one parallel 
to the diagonal and 



\\R{z;s)\\=sup{\R(z) i% 



3\> 



J 



s} < c2-\ s y 



Hence the series (139|) converges in the operator norm. With some abuse of notation, 
we shall denote the corresponding bounded operator again by the symbol R(z). 
Further one observes that, on the level of formal matrix products, 



Q-z)R{z)=R{z){3-z) = I. 
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The both equalities are in fact equivalent to the countable system of equations (with 
w = 0) 



w k -\G(z; i, h — 1) + (A fc - z)G(z; i, k) + w k G(z; i, k + 1) = 5 itk , i, k G N. 

This can be verified, in a straightforward manner, with the aid of the rule (j2J) or some 
of its particular cases (j^J) and (jSJ). By inspection of the domains one then readily 
shows that the operators J max — z and R(z) are mutually inverse and so z G p(Jmax)- 
Finally, suppose that |A n | — > oo as n — > oo, and z G C \ 3(c7)- It turns out that 
then the above estimates may be somewhat refined. In particular, (|38|) is majorized 

by 

| Ai -*|-V2| A ._ z |-l/2 2-i-M 

for n < But this implies that R(z; s)jj — >■ as z, j — >• oo, with i — j = s being 
constant. It follows that the operators R(z; s) are compact. Since the series 
converges in the operator norm, R(z) is compact as well. □ 

Corollary 13. If condition (TjSp is fulfilled for at least one z G Cq then 

spec(J max ) \ der(A) = spec p (J max ) \ der(A) = 3(J)- 
Proof. Propositions [TT] and [T2] respectively imply the inclusions 

3(J r ) C spec p (J max ) \ der(A), spec(J max ) \ der(A) C 3(J)- 
This shows the equality. □ 



Theorem 14. Suppose that the convergence condition [28)) is fulfilled for at least one 
z G Cq and the function Fj(z) does not vanish identically on Cg. Then J min = 

spec(J) \ der(A) = spec p (J) \ der(A) = 3(J). (40) 

Suppose, in addition, that the setC\ der(A) is connected. Then spec( J)\der(A) consists 
of simple eigenvalues which have no accumulation points in C \ der(A). 

Proof. By the assumptions, C \ (der(A) U 3(i7")) 7^ 0- From Proposition [T21 one infers 
that p(J max ) 7^ 0- According to Theorem [101 one has J m ; n = J max . Then (j4"0|) becomes 
a particular case of Corollary [13j 

Let us assume that C\ der(A) is connected. Then the set Cq is clearly connected 
as well. Suppose on contrary that the point spectrum of J has an accumulation point 
in C \ der(A). Then, by equality (1401) . the set of zeros of the analytic function Fj(z) 
has an accumulation point in C \ der(A). This accumulation point may happen to be 
a member A„ of the sequence A, but then one knows that Fj(z) has a pole of finite 
order at A n . In any case, taking into account that Cq is connected one comes to the 
conclusion that Fj(z) = everywhere on Cq, a contradiction. □ 
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Remark 15. Theorem [TH is derived under two assumptions: 

(i) The convergence condition (|28|) is fulfilled for at least one Zq G Cq. 

(ii) The function Fj(z) does not vanish identically on Cq. 

But let us point out that assumption (ii) is automatically fulfilled if (i) is true and 
the range of the sequence A is contained in a halfplane. This happens, for example, 
if the sequence A is real or the sequence {ReA n }^ =1 is semibounded. In fact, let us 
for definiteness consider the latter case and suppose that Re A n > c, Vn G N. Then 
(— oo, c) C Cg and l/|A n — z\ tends to monotonically for all n as z — > — oo. Similarly 
as in (131) one derives the estimate 



|iM*)-l| <exp [J2 



v n=l 



w 2 



(A„ — z)(X n+ i — z) 



- 1. 



It follows that \im z ^_ OQ Fj(z) = 1. Notice that in the real case, the function Fj(z) 
can identically vanish neither on the upper nor on the lower halfplane. 

Corollary 16. Let J be real and suppose that [23\) is fulfilled for at least one Zq G Cq. 
Then J min = J max = J is self-adjoint and spec( J) \ der(A) = 3(i7) consists of simple 
real eigenvalues which have no accumulation points in M. \ der(A). 

Proof. Some assumptions in Theorem [TH become superfluous if J is real. As observed 
in Remark [TBI assuming the convergence condition the function Fj{z) cannot vanish 
identically on Cq. The operator J is self-adjoint and may have only real eigenvalues. 
The set C \ der(A) may happen to be disconnected only if the range of the sequence A 
is dense in M, i.e. der(A) = R. But even then the conclusion of the theorem remains 
trivially true. □ 

Let us complete this analysis by a formula for the norms of the eigenvectors de- 
scribed in Proposition [TTJ In order to simplify the discussion we restrict ourselves to 
the domain Cq. Then instead of fl32|) one may write 



This is in fact nothing but the solution f n from (121|) restricted to nonnegative indices. 
Proposition 17. If z G Cg satisfies then the functions £,k( z )> k G Z + , defined in 

B fulfill 

oo 

E^W 2 = ^(^W-^(^W- ( 42 ) 
fe=i 

Particularly, if in addition J is real and z G R R Cq is an eigenvalue of J then 
£( z ) = (£k( z ))kLi i s a corresponding eigenvector and 

U(z)\\ 2 = Z'o(zMz)- (43) 
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Proof. Put (j 1 ^ = z — Xj, (j 2 ^ = y — \j, j G N, in equation (123|) . where Cg. Then 

Proposition [6] is applicable to /j 1 ^ = £j(z), ff^ = £j(y), j G Z+. Hence (u>o = 1) 

oo 

(2 - y) £>(*)&(!/) = &(*)6>(j/) - £o(*)6(y)- 

fc=0 

Now the limit y z can be treated in a routine way. □ 

Corollary 18. Suppose J is real and let condition |15J) be fulfilled for at least one 
z o £ Q). TTien £/ie function Fj(z) has only simple real zeros on 

Proof. Suppose z G Cg is a zero of i.e. Fj-(z) = £o(z) = 0. Then z is a real 

eigenvalue of J where J = J max = J m j n is self-adjoint, as we know from Corollary H"6l 
Moreover, by Proposition [TTJ 7^ is a corresponding real eigenvector. Hence from 
(143 p one infers that necessarily ^l- 2 ) 7^ 0- ^ 



3.4 Approximation of the spectrum by spectra of truncated 
matrices 

Let us first introduce some additional notation which will be needed in the current 
subsection. W and L stand for the diagonal matrix operators whose diagonals equal 
w and A, respectively. U designates the unilateral shift and U* its adjoint operator 
(Ue n = e n+ \ for n = 1, 2, 3, . . . , with e n being the vectors of the standard basis). 
For a Jacobi matrix J we introduce the set 

A(J) ■= e C; lim dist(spec(J n ), //) = 0} (44) 

n— >oo 

where J n is defined in fTTOT) but now with v = w. Thus fi G A(j7") iff there exists a 
sequence of eigenvalues {/!„,} of J n such that jj, n — > \i as n — >■ 00. 

Remark 19. (i) Definition (|44|) is taken, for example, from [3]. Some authors prefer 
to work, however, with the set A(j7") formed by all limit points of sequences {/!&} of 
eigenvalues from spec(J nk ), with {n^} being any possible strictly increasing sequence 

of indices [TT | IT2" ] . One clearly has A(J~) C A(J') and the inclusion is in general strict 

as demonstrated by a simple example in [TT| Eq. (2)]. 

(ii) In the case where the sequences A and w are real and positive, respectively, 
the relation between the spectrum of a Jacobi operator J and the set A( c 7) has been 
intensively studied by several authors. Most results of this kind are restricted to the 
bounded case, however. The inclusion spec(J) C A(jT') is proved for bounded Jacobi 
matrices in [3J, Theorem 2.3]. It cannot be replaced, in general, by the sign of equality 
as demonstrates an counterexample constructed in the Appendix in [3J. 

(iii) In [T4"| Theorem 5.1] it is shown that spec(Js) C A(J~) where the operator 

J S :=L + WU* + UW 

(an operator sum) is assumed to be self-adjoint. Some reasoning used below, notably 
that in Lemmas |2"T1 and I2"2"j is inspired by this article. We also remark that Theorem |2"B1 
below partially reproduces and overlaps with Theorems 2.1 and 2.10 from [12J. 



19 



Lemma 20. If condition / f^) is fulfilled for at least one Zo G Cq then 

A(J)\dex(X) c3(J). 



(45) 



Proof. Let z G A(j7") \ der(A). Then, by definition, there exists a sequence {fJ, n } such 
that fj, n G spec(J„) and /i n — > z as n — > oo. Considering sufficiently large indices n 
one may assume that fi n does not coincide with any member Afc of the sequence A, 
possibly except of z. Using (fT3|) one gets, for all sufficiently large n, 



( 



= det(fi n I n - J n ) 



\ 



\\ {f^n - Afc) 



\ k=l 



/ 



lim (u — z 



r(z) 



Xi-u 



1=1 



(46) 



Denote temporarily 

F n {u) = {u-z) r ^$ 



7/ 



Xi-u 



i=i 



, u G Cq, JiGNU {oo}. 



All functions F n (u) have a removable singularity at u = z. Moreover, with some slight 
modification of Lemma M one can show that F n (u) — > F°°(u) as n — > oo uniformly 
on a neighborhood of the point z. Since all the involved functions are analytic on 
this domain, one knows that the corresponding derivatives converge uniformly as well. 
Equation (146]) means that F n (/i n ) = for all n sufficiently large. One concludes that 



lim (u - z) r{z) Fj(u) = F°°(z) = lim F"(/i T 



0. 



Hence z G 3(c7)- 



□ 



Denote by P n , n G N, the orthogonal projection on £ 2 (N) onto the linear hull of 
the first n vectors of the standard basis. Then the finite-dimensional operator J n 
introduced in (1101) can be identified with P n JP n restricted to the subspace RanP n . 
Here J is any operator such that J C J (see Definition 

Lemma 21. Let J be any operator on £ 2 (N) fulfilling J C J. Then (in £ 2 (N)) 
lim P n JP n x = Jx, Vx G Dom(J) n Dom(W*). 

Proof. Let x G Dom(J) H Dom(WU*). Since 

(J ~ P n JPn)x = (I- P n )Jx + P n J(I - P n )x 

one has 



\\(J-P n JP n )x\ 



(I-P n )Jx\\ 2 + \w n \ 2 \(x,e n+1 )\ 2 = \\(I-P n )Jx\\ 2 + \(WU*x,e n ) 



The RHS clearly tends to zero as n — > oo. 



□ 
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Lemma 22. Suppose that J is real and let J be a self-adjoint extension of the sym- 
metric operator J m in- //Dom(J) C Dom(WU*) thenspec(J) C A(J). If, in addition, 
condition |J2|) is fulfilled for at least one zq G Cq then 

3G7) = spec(J) \ der(A) = spec p (J) \ der(A) = A{J) \ der(A). (47) 

Proof. First suppose that, on the contrary, fi G spec(J) and fi (ji A(J). Then there 
exist d > and a subsequence {J nk }^ =l such that dist(//, spec J nk ) > d for all k. Let 
x be an arbitrary vector from Dom(J) C Dom(WU*). 

As is well known and easy to verify, if A is a self-adjoint operator such that 
belongs to its resolvent set then ||^4 _1 || < 1/ dist(0, spec A). This implies that 

V/GDomi, \\Af\\ >dist(0,specA)||/||. 

Applying this observation to the operators J nk one gets 

\\{fiPn h - Pn k JPn k )x\\ > d\\P nh x\\. 

Sending k — > oo and referring to Lemma I2T1 one concludes that — J)x\\ > 
for every x G Dom(J). By the Weyl criterion fi G p{J), a contradiction. 

If, in addition, (j28|) is fulfilled for at least one Zq G Cq then we know from Corol- 
lary dn] that J is in fact unique and the first two equalities in (1471) are true. So 
combination of Lemma [2U] with the first part of the current lemma yields 

3(J) = spec(J) \ der(A) C A(J) \ der(A) C 3(J). 

This completes the proof. □ 

Theorem 23. Suppose that the Jacobi matrix J is real and let condition (d5j) be 
fulfilled for at least one z 6 Cq. If any of the following two conditions is satisfied: 

(i) XmVn^Wn = 0, 

(ii) Hindoo |A n | = oo and 

limsup j— + limsup p- —. < 1, (48) 

n— >oo \^n\ n—>oo l^n+ll 

then 

A(J) = spec(J), (49) 

where J is the unique self-adjoint operator whose matrix in the standard basis equals 
J. 

Proof. According to Corollary [TSJ J = J min = J max is the unique self-adjoint operator 
determined by J . 

(i) Let lircin^ooWn = 0. Then we claim that 

der(A) = spec ess (J) C spec(J) C A(J). (50) 
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In fact, the Hermitian operator WU* + UW is compact and 

J — L + (WU* + UW) (51) 

(an operator sum, Dom J = DomL). But the Weyl theorem tells us that J has the 
same essential spectrum as L, which is nothing but der(A). The last inclusion in (150]) 
follows from Lemma 1221 Again by Lemma 122] spec(J) \der(A) = A(j7") \der(A), which 
jointly with (150]) implies (H9|) . 

(ii) Suppose that |A„| — > oo and f ]4"5]) is true. Then clearly der(A) = 0. We wish to 
verify the assumptions of Lemma 1221 Since the sum of a self-adjoint operator with a 
bounded Hermitian operator does not change the domain one may suppose, without 
loss of generality, that 

i I \ i w i \ W n\ \w n \ ^ 

I A n | > 1, Vn, and sup — — - + sup — < a < 1, 

n \^n\ n |^n+l| 

where a is a constant. But this means that 



|| (W + UW)L~ l \\ < \\WU*L~ 1 \\ + \\WL^\\ < 



and hence (WU* + UW) is L-bounded with a relative bound smaller than 1. By the 
Kato-Rellich theorem [TTj Theorem V.4.3], (15T|) is again true. Moreover, 



Dom J = DomL C Dom(W*). 

Hence in this case, too, Lemma [221 implies (149]) . 
Remark 24. Observe that (|2"8|) already implies that 



□ 



lim 

n— ¥oo 



w; 



0. 



Thus assumption 



A n A n +i 

may well turn out to be superfluous in some concrete examples. 



4 Examples 

4.1 Explicitly solvable examples of point spectra 

In all examples presented below the Jacobi matrix J is real and symmetric. The set 
of accumulation points der(A) is either empty or the one-point set {0}. Moreover, 
condition f |28|) is readily checked to be satisfied for any zq e C \ WL. Thus Corollary [16] 
applies to all these examples and may be used to determine the spectrum of the unique 
self-adjoint operator J determined by J (recall also definition f]3T]) of the zero set of 
the characteristic function). In addition, Proposition [TT1 and equation ( 152"]) (or (|4l|) ) 
provide us with explicit formulas for the corresponding eigenvectors. 
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Example 25. This is an example of an unbounded Jacobi operator. Let X n = na, 
where a£l \ {0}, and w n = w > for all n G N. Thus 



J 



(aw \ 

w 2a w 

w 3a w 

V '-J 



One has der(A) = and spec(J) = ~5{J). Using f fT5"j) one derives that 



5 



ak — z 



k=r+l 

for r G Z_i_. It follows that 



a/c — 2 



fc=r+l 



rfl + r- -]J r _ z/a ( —J 



a. 



spec(J) = |z6C; J_ 2 / Q ( —J = 



2w\ 



(52) 



For the corresponding eigenvectors v (z) one obtains 

Vk ( z ) = (-i) k j k _ z/a (^y kem 

Let us remark that the characterization of the spectrum of J, as given in 
was observed earlier by several authors, see |15[ Sec. 3] and |19[ Thm. 3.1]. We 
discuss in more detail solutions of the characteristic equation J_ 2 (2w) = below in 
Subsection 14.31 

Further we describe four examples in which the Jacobi matrix always represents a 
compact operator on £ 2 (N). We shall make use of the following construction. Let us 
fix positive constants c, a and (3. For n G Z + we define the c-deformed number n as 



n 



n-l 
i=0 



Hence [n] c = (c n — l)/(c — 1) if c ^ 1 and [n] c = n for c = 1. Notice that 
[n + m — l] c — [n — l] c 



He 

As for the Jacobi matrix we put 

1 



[n] c — [n — l] c , Vn, m G N. 



(53) 



A, 



a; + [n - l] c 



w 



n — P\/^n ~ A n+ i , n — 1,2,3,.... 



(54) 
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We claim that the identity 

oo oo oo 

E E - E 

k±=r fc2=fci+2 fc s =fc s _i+2 



X 



W 



k 2 



"I)' 



(A fel - z)(\ kl+1 - z) (X k2 - z)(\ k2+1 - z) (X ks - z)(X ks+1 - z) 
2S fl([^-T^-)Y (55) 



i=l 



holds for every r,s G N. In fact, to show (|55|) one can proceed by mathematical 
induction in s. The case s = 1 as well as all induction steps are straightforward 
consequences of the equality 



w i 



P 2 



n I* 1 



-i 



n [-Mi 



with s = 1, 2, 3, . . ., which in turn can be verified with the aid of (l53]h 
A principal consequence one can draw from ( 155]) is that 

,2 1°° \ 00 o2s 



3 



in 



E^n([>M> 



An Z ) n=r+1 J s=Q Z* r t 



A 



r+i 



(56) 



holds for all r G Z+. 

Example 26. In (I54p . let us put c = 1 and a = 1 while (3 is arbitrary positive. Then 
A„ = 1/n, W7„ = /?/ ^/n(n + 1) , for all nGl, and so 

/ 1 0/V2 \ 
P/y/2 1/2 P/VE 

P/y/E 1/3 /3/\/l2 



J 



V 



One finds that 

2 oo 







0^(^ + 1--;- — 



(57) 



with G Z+, see UJ Eq. 9.1.69]. Then 



*w = n i - 1) (f 
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-l/a 



/ 1 2/3 



and 



spec(J) ={zeR\ {0}; J_ 1/jB ( ^ j = } U {0}. 



For the corresponding eigenvectors t> (2;) one has 



Example 27. Now we suppose in ()54|1 that c > 1 and put a = l/(c — 1). Then 
A n = (c— l)c _n+1 and w n = (3 (c— l)c -1 / 2 c(~ n+1 " 2 . In order to simplify the expressions 
let us divide all matrix elements by the term c — 1. Furthermore, we also replace the 
parameter (3 by /3c 1//2 and use the substitution c = 1/g, with < q < 1. Thus for the 
matrix simplified in this way we have A n = q n ~ l and w n = /3g( n_1 " 2 . Hence 



J 



/l /3 \ 

/3 q (3^q 

(3^/q q 2 0q 

v '•■ '•■ '■) 



Using (I17p . equation ( 15 6 p then becomes 

2 -\ 00 



/ Ti 



X n — z 



n=r+l 



s q s ^ 



{q;q)s(q r /z;q) s \ z 



2.s 



q r /z;q,-q r /3 2 /z 2 ), r G Z + 



Thus we get 



spec(J) = j* £l \ {0}; ; g 



o<Pi[\^q,~) = 0|>U{0}. 



The fcth entry of an eigenvector v(z) corresponding to a nonzero point of the spectrum 
z may be written in the form 

V k (z) = g(*-D(*-^ f ^V _1 ; g) fc( ; £ ; g, ) , k G N. 



z / 00 



Further we shortly discuss two examples of Jacobi matrices with zero diagonal 
and w G ^ 2 (N). Such a Jacobi matrix represents a compact operator (even Hilbert- 
Schmidt). The characteristic function is an even function, 



<{?}") -e^Pe e 

\K. ) n=i/ m=0 fel=1 fca=fcl . 



E 



+2 fc m — fc m _i+2 



Hence the spectrum of J is symmetric with respect to the origin. 

Though always belongs to the spectrum of a compact Jacobi operator, one may 
ask under which conditions is even an eigenvalue (necessarily simple). An answer 
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can be deduced directly from the eigenvalue equation fTTTh . One immediately finds 
that any eigenvector x must satisfy x 2 fc = and %2k-i = (— l) fe+1 2 ; i/72fc-i, k G N. 
Consequently, zero is a simple eigenvalue of J iff 

fc=1 72*-l 



Example 28. Let A n = and w n — 1/ y^(n + a)(n + a + 1), n£N, where a > — 1 is 
fixed. According to (ITS]) one has, for k G Z + , 

<-y 2 1 °° \ / f 1 1 °° \ /2 



Hence 

spec(J) = jz G M\{0}; J a Qj = J U {0}. 

The fcth entry of an eigenvector v(z) corresponding to a nonzero eigenvalue z may be 
written in the form 



(z) = Va + k z a J a+k 0J , fcGN. 



It is well known that for a G 1/2 + Z, the Bessel function Jq,(z) can be expressed 
as a linear combination of sine and cosine functions, the simplest cases being 

J-i/2(z) = J— cos(z), Ji/ 2 (z) = \ — sin(z). 

V VT2 V 

Thus for a = ±1/2 the spectrum of J is described fully explicitly. In other cases 
the eigenvalues of J close to zero can approximately be determined from the known 
asymptotic formulas for large zeros of Bessel functions, see (U Eq. 9.5.12]. 

Example 29. Suppose that < q < 1 and put A n = 0, w n = q n ~ l , n G N. With the 
aid of (11 7p one derives that 

It follows that 

spec(J) = {zeR\ {0}; 0; g 2 , -^T 2 ) = 0} U {0}. 

The components of an eigenvector v(z) corresponding to an eigenvalue z ^ may be 
expressed as 

Vk(z) = gC- 1 ^^ ^( ; 0; g 2 , -g 2fc ^ 2 ), fc G N. 

In this example as well as in the previous one, belongs to the continuous spectrum 
of J since the condition (158)) is not fulfilled. 
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4.2 Applications of Proposition [171 

Here we apply identity fj42l) to the five examples of Jacobi matrices described above. 
Without going into details, the final form of the presented identities is achieved after 
some simple substitutions. On the other hand, no attempt is made here to optimize 
the range of involved parameters; it is basically the same as it was for the Jacobi 
matrix in question. 

1) In case of Example 1251 one gets 

fc=l ^ 

where x > and i/gC. This is in fact a particular case of f llSp . 

2) In case of Example 1261 one gets 

°° z 2 ( d d \ 

k J- az +k{z) 2 = — ( J- az (z) — J- az+ \{z) — J-az+l{z) J^^ 02 ^) ) 
k=l \ Z Z J 

where a > and zGC. 

3) In case of Example |2"51 one gets 

°° z 2 ( d d \ 

^2 ( a + k)J a+k (z) 2 = — I J a (z) — J a+l (z) - Ja+l{z) ^ J a( Z ) J > 
k=l ^ ' 

where a > — 1 and z £ C. 

4) In case of Example |27] one gets 

oo 

ST q ^-m-m {(3z) ^ (tf z . qU qM . gk z . Qj - q K(3* z 2 )) 2 



k=l 



(Q z ; ( O0i(; Z] q, -fi 2 z 2 ) o 0i(; qz; q, -qf3 2 z 2 ) 

( d 
+ *(*-!) (o<M; g^; g, -q(3 2 z 2 ) — <M; g, -/?V) 



- o<M; «; g, ^ o<M; g^; g, -g/?Vm 



where 0<g<l,/3>0 and zeC. 
5) In case of Example 1291 one gets 



oo 

j^^-^^^o^iOO;?,-?**) 8 = o0i(;O;g,-^) o 0i(;O;g,-g^) 

k=l 

+ 2z^ o 0i(;O;g, -2) ^ o 0i(;O;g, -gz) - <M; 0; g, ^ o 0i(; 0; g, 



where < q < 1 and z G 
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4.3 A Jacobi matrix with a linear diagonal and constant par- 
allels 



Here we discuss in somewhat more detail Example [25] concerned with a Jacobi matrix 
having a linear diagonal and constant parallels. For simplicity and with no loss of 
generality we put a = 1. Our goal is to study how the spectrum of the Jacobi 
operator J depends on the real parameter w. We treat J as a linear operator- valued 
function, J = J(w). One may write J(w) = L + wT where L is the diagonal operator 
with the diagonal sequence A n = n, Wi G N, and T has all units on the parallels 
neighboring to the diagonal and all zeros elsewhere. Notice that ||T|| < 2. 

We know that J(w) has, for all w G R, a semibounded simple discrete spectrum. 
Let us enumerate the eigenvalues in ascending order as X s (w), s G N. From the 
standard perturbation theory one infers that all functions X s (w) are real analytic, 
with A s (0) = s. Moreover, the functions X s (w) are also known to be even and so we 
restrict w to the positive real half-axis. In Example [25] we learned that for every w > 
fixed, the roots of the equation J^ z (2w) = are exactly X s (w), s G N. Several first 
eigenvalues X s (w) as functions of w are depicted in Figure 1. 

The problem of roots of a Bessel function depending on the order, with the ar- 
gument being fixed, has a long history. Here we make use of some results derived in 
the classical paper [6]. Some numerical aspects of the problem are discussed in |15j . 
For comparatively recent results in this domain one may consult [19] and references 
therein. 

In [6J it is shown that 



From this relation one immediately deduces a few basic qualitative properties of the 
spectrum of the Jacobi operator. 

Proposition 30 (M. J. Coulomb). The spectrum {X s (w); s G N} of the above intro- 
duced Jacobi operator, depending on the parameter w > 0, has the following properties. 

(i) For every s G N, the function X s (w) is strictly decreasing. 

(ii) If r < s then X r '(w) < X s '(w). 

(Hi) In particular, the distance between two neighboring eigenvalues X s+ \(w) — 
X s (w), s 6 N, increases with increasing w and is always greater than or equal to 1, 
with the equality only for w = . 

Let us next check the asymptotic behavior of X s (w) at infinity. The asymptotic 
expansion at infinity of the sth root j s (y) of the equation J v {x) = 0, with v being 



where a s is the sth negative zero of the Airy function Ai(x). From here one deduces 
that 




fixed, reads [H Eq. 9.5.22] 



Jsiy) =v-2 



/3 a s v 1/3 + 0(v 




as v — > +oo, 




as w — > +oo. 
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Concerning the asymptotic behavior of X s (w) at w — 0, one may use the expression 
for the Bessel function as a power series and apply the substitution, X s (w) = s — z(w), 
s = 1, 2, 3, . . .. The solution z = z(w), with z(0) = 0, is then defined implicitly near 
w = by the equation 

Y ^ w 2m = 0. 

^— ' ml T(m + 1 — s + z(w)) 

771=0 

The computation is straightforward and based on the relation 

r , l — ; = ("l) m m! (z - ^ (0) (m + l)z 2 ) + 0(^ 3 ) , m = 0, 1, 2, 3, . . . , 
1 (— m + zj v / \ / 

where ip^ is the polygamma function. This way one derives that, as w — > 0, 
X^w) = l- w 2 + -w i + 0(w 6 ), (59) 

X s (w) = S-- 1 -— W 2s + - -A- w 2 S+2 + / 2S+ 4X fors>2 . 

The same asymptotic formulas, as given in ( 15 9 p . can also be derived using the 
standard perturbation theory |T7l § II. 2]. Alternatively, one may use equivalent for- 
mulas for coefficients of the perturbation series derived in [TJ [8] which are perhaps 
more convenient for this particular example. 

The distance of s G N to the rest of the spectrum of of the diagonal operator L 
equals 1. The Kato-Rellich theorem tells us that there exists exactly one eigenvalue 
of J(w) in the disk centered at s and with radius 1/2 as long as \w\ < 1/4. The 
explicit expression for the leading term in (159]) suggests, however, that the eigenvalue 
X s {w) may stay close to s on a much larger interval at least for high orders s. It 
turns out that actually X s (w) is well approximated by this leading asymptotic term 
on an interval [0, /3 a ), with f3 s ~ s/e for s 3> 1. A precise formulation is given in 
Proposition [23 below. 

Denote by yk{v) the /cth root of the Bessel function Y u (z), fceN. Let us put 

' s -l)!^ 1/(2s) 



0s ■= ^~ ^—J , s G N. 

In order to avoid confusion with the usual notation for Bessel functions, the nth 
truncation of J(w) is now denoted by a bold letter as J n (w). 

Lemma 31. The following estimate holds true: 

Ps<lyi(s-\\ Vsgn. (60) 

Proof. One knows that v < yi(v), W > (U Eq. 9.5.2], and in particular this is true 
for v = s — 1/2, s G N. On the other hand, the sequence 



l)!s! 
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is readily verified to be increasing, and 1 < 04. This shows (160|) for all s > 4. The 
cases s = 1,2,3 may be checked numerically. □ 

Lemma 32. Denote by Xn{w; z) the characteristic polynomial of the nth truncation 
J n (w) of the Jacobi matrix J(w). IfO < w < (3 S for some seN then z = X s (w) solves 
the equation 

wJ 2 s~z(2w) 

X2s-i(w; z) - —— X2s _ 2 (w;z) = 0. (61) 

J 2 s-l-z{2w) 

Proof. Let {e^; k G N} be the standard basis in £ 2 (N). Let us split the Hilbert space 
into the orthogonal sum 



£ 2 {N) = span {e k ; 1 < k < 2s - 1} © span {e k ; 2s < k}. 
Then J(w) splits correspondingly into four matrix blocks, 

MH B( W ) 

J{W> \ C(w) D(w) 

Here A(w) = J2 S -i{w), D(w) = J{w) + (2s — 1)/, the block B(w) has just one nonzero 
element in the lower left corner and C(w) is transposed to B(w). 

By the min max principle, the minimal eigenvalue of D(w) is greater than or equal 
to 2s — 2w. Since X s (w) < s one can estimate 

min spec (D(w)) — \ s (w) = \i(w) — \ s (w) + 2s — 1 > s — 2w. 

We claim that < w < j3 s implies minspec(D(w)) — X s (w) > 0. This is obvious for 
s = 1. For s > 2, it suffices to show that (3 S < s/2. This can be readily done by 
induction in s. Hence, under this assumption, D(w) — z is invertible for z = X s (w). 

Solving the eigenvalue equation J(w)v = zv one can write the eigenvector as a 
sum v = x + y, in accordance with the above orthogonal decomposition. If D(w) — z is 
invertible then the eigenvalue equation reduces to the finite-dimensional linear system 

(A- z- B(D - z)- 1 C)x = 0. (62) 

One observes that B(D — z)~ 1 C has all entries equal to zero except of the element 
in the lower right corner. Using (135|) and (fT5l) one finds that this nonzero entry equals 

wJ 2 s~z(2w)J 2 s-i-z(2wy 1 . 

Equation (16T|) then immediately follows from (162]) . □ 

Proposition 33. For s E N and < w < (3 S , one has 

\ / \ 1 / mv 2s 

< s - X s (w) < - arcsin -— - 

7T V s — 1 !s! 
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Proof. We start from Lemma [32] and equation f }6T|) . Let us recall from [20, Proposi- 
tion 30] that 



det (J 2s -i(w) -a-x) = (-l) s x ^ 
Hence if 2: £ K, \z — s\ <1, then 



S X / r-\ 7 -t \ S n> 1 

2s — k — 1 



k=0 



w 2k "Q 

3=1 



- „2j 



s-l 



Ix2*-i (w; «)| > \z - s\ Yl (f - - s ) 2 )- 

3 = 1 



(63) 



Since J_ s+1 / 2 (:r) = (— l) s 5^-i/2(^) it is true that for 2w = yi(s — 1/2) one has 
X s (w) = s — 1/2. Because of monotonicity of X s (w) one makes the following observa- 
tion: if 2w < yx{s - 1/2) then s > X 3 (w) > s - 1/2. 

By Lemma l3T| if w < /3 a then 2to < — 1/2), and so the estimate (163]) applies 
for z = A s (w). Using also Proposition 2] to express x 2s „ 2 (w;; z) one derives from flBT]) 
that 



|A — s| < w 



J 2s _a(2w) 



J 2s -i-a(2w) 



s-X 



w 



w w 

2s - 1 - A u V l^A' 2 - A' " ' ' 2s - 2 - A 



(64) 



where as well as in the remainder of the proof we write for short A instead of X s (w). 
Starting from the equation 



w w 



w 



1-A'2-A'3-A' 
and using (jlj, f[T5"]) one derives that, for all k G Z + , 



0, with A = XJw) 



WW w 

? ■ ■ ■ J 



A'2-A' 'jfe-A 



fc J fc+ i- A (2w) 
Ji-a(2w) 



(65) 



Combining ( |64"]) and ( 165]) we get (knowing that < s — A < 1/2 for A = A s (u>)) 



s — A < w 



2s-l 



-1 s-l 
V?°=l 2=1 



J 2s _ A (2w) 



Ji-a(2w) 

But notice that, by expressing the sine function as an infinite product 



s-l 



3=1 

Hence 



3=1 



n 



sin(7r(s — A)) < n 



sin(7r(s — A)) 
7r(s - A) 

w 28 - 1 J 2s -x(2w) 



\3=s 



J 2 



((s-l)O 2 
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Ji-x(2w) 



From (126|) one gets, while taking into account that J-\(2w) = 0, 

sin(7rA) = ttwJ\(2w) J\-\(2w). 
In addition, one knows that 

1 



\Ju x ) < 



I> + 1 



X 



provided v > -1/2 and s 6 1 jH Eq. 9.1.62]. Hence 
sin(7r(s - A)) 2 < 



71 2 W 4S 



((s-l)!) 2 r(2s + l-A)r(A + l) ' 
Writing A = s — £, with < ( < 1/2, one has 

1 



-1 

Thus we arrive at the estimate 



^ w (s + C + i) + # 0) (s-C + i) <o. 



sin(7r(s - A)) 2 < 



7T 2 W 4s 



((,-l)!)2( S !) 2 



To complete the proof it suffices to notice that the assumption w < j3 s means 
nothing but w 2s / ((s — l)!s!) < 1, and it also implies that < s — A < 1/2. □ 
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Figure 1. Several first eigenvalues X s (w) as functions of the parameter w for the Jacobi 
operator J = J(w) from Example [25] with a — 1. 
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